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During  the  fourth  year  of  the  contract,  further  advances  were  Bade 
towards  the  goal  of  gathering  the  heat  transfer  and  aerodynamics  flow  data 
necessary  for  a good  understanding  of  the  performance  of  film-cooled, 
highly-loaded,  transonic  turbine  blading. 

Surface  Mach  number  and  heat  transfer  rate  distributions  were 
determined  for  a reference  transonic  airfoil  over  a range  of  exit  Mach 
numbers,  for  inlet  incidence  angle  variation  of  ±15°.  An  evaluation  and 
comparison  of  all  cascade  data  collected  so  far  was  then  conducted. 
Progress  was  also  made  in  the  investigation  of  the  effects  of  unsteadiness 
on  transonic  airfoil  aerodynamics  and  heat  transfer. 
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STUDIES  ON  TRANSONIC  TURBINES  WITH 


FILM-COOLED  BLADES 
A FOURTH  ANNUAL  REPORT 

N.  Adana,  F.  Ha j jar,  R.  F.  Topping,  J.  F.  Louia 

Introduction 

i 

i 

The  increaaing  intereat  in  tranaonic  turbinea  ahown  over  the  laat  few  , 

l 

years  is  a reflection  of  the  need  to  increase  turbine  work  output  per  unit 
area  of  cooled-blade  surface.  Much  of  the  increase  in  turbine  cycle 
efficiency  obtained  over  the  laat  two  decades  has  been  due  to  increased 
turbine  inlet  temperature  through  the  use  of  more  efficient  cooling 
techniques.  As  more  and  more  coolant  mass  flow  is  extracted  from  the 
compressor,  however,  a point  is  reached  where  the  losses  and  reduction  in 
working  fluid  caused  by  the  use  of  increased  coolant  mass  flow  more  than 
offset  the  increase  in  efficiency  and  specific  thrust  attainable  by  higher 
turbine  inlet  temperatures . At  this  point,  the  need  arises  to  increase  and 
optimize  the  work  output  per  unit  of  cooled-blade  surface.  This  can  be 
accomplished  by  the  use  of  transonic  blading. 

The  goal  of  these  studies  is  to  gather  necessary  heat  transfer  and 
aerodynamic  data  for  the  designer  of  high  performance,  film-cooled, 
highly-loaded  transonic  turbine  blading.  In  the  first  three  years  of  the 
contraot,  the  following  key  tasks  were  performed: 
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(a)  The  gathering  of  film  cooling  effectiveness  data  in  the 
transonic  range  1.15<M<1.4  for  both  slot  and  hole 
injection  using  a shock  tunnel. 

(b)  The  gathering  of  heat  transfer  and  aerodynamics  data  at 
the  trailing  edge  of  a transonic  blade,  and  the  use  of 
this  data  to  formulate  a model  for  the  heat  transfer 
and  aerodynamic  flow  at  the  trailing  edge  of  transonic 
blades. 

i 

(e)  The  design,  manufacture  and  check-out  of  the  MIT 
cascade  blowdown  facility. 

(d)  The  design  and  aerodynamic  testing  of  four  transonic 
blade  profiles  in  a conventional  wind  tunnel  at  VKI. 

(e)  Aerodynamic  and  heat  transfer  testing  of  the  four  I 

transonic  blade  profiles  at  zero  angle  of  incidence  in 

the  blowdown  facility. 

(f)  The  design  and  manufacture  of  equipment  and 
instrunentation  to  measure  the  effects  of  periodic 
unsteadiness  on  the  aerodynamic  and  heat  transfer 
performance  of  turbine  blades. 

(g)  An  estimation  of  the  coolant  flow  requirements  of  a 
transonic  turbine  and  a comparison  between  the  coolant 
flow  requirements  of  the  transonic  turbine  with  that  of 
a subsonic  turbine  of  equal  work  output. 


-2 


In  the  fourth  year,  the  following  tasks  were  accomplished: 

(1)  Aerodynamic  and  heat  transfer  tests  were  performed  on 
the  reference  blade  profile  at  off-design  angles  of 
incidence  in  the  blowdown  facility. 

(2)  Heat  transfer  and  aerodynamic  properties  of  the 
reference  blade  were  compared  with  earlier  experimental 
work. 

(3)  The  experimental  and  analytical  program  on  the  effects 
of  unsteadiness  on  transonic  airfoils  began. 


r 

■ 

1.  Aerodynamic  and  Heat  Transfer  Testing  at  Off -Design  Incidence  ] 

1-1  The  Effect  of  Incidence  on  the  Pressure  Distribution  Around  The j 

j Reference  Blade  j 

j ;! 

j Figure  1 shows  the  reference  blade  used  in  conducting  the  incidence  ] 

i I 

tests  and  shows  the  pertinent  parameters  of  the  blade.  Table  1 shows  the  \ 

■ } 

location  of  the  pressure  taps  on  the  blade,  the  distance  x being  measured  j 

,,  ' i 

J.  -J 

i along  the  surface  of  the  blade.  The  isentropic  surface  Mach  numbers  were  ] 

( ) 

' plotted  as  a function  of  the  non-dimensional  length  ratio,  S/L,  for  the 

\ I 

[■  design  incidence  of  0°  and  off-design  incidences  of  ± 7.5  and  ±15°.  The 

* range  of  the  isentropic  exit  Mach  number  varied  from  the  subsonic  region  ■ 

[ through  the  low  supersonic  region.  The  reference  blade  was  the  same  as 

that  used  in  previous  work  with  the  exception  that  surface  roughness  was 

l 

r V 

( removed  by  smoothing  the  4X  fabrication  master  blade. 


Figures  2a  through  2e  are  plots  of  the  Mach  number  distribution  around 
the  blade  surface  for  increasing  mass  flows.  The  stagnation  point  was 
found  using  a potential  flow  analysis.  A line  perpendicular  to  the  flow 
angle  was  used  to  find  the  tangent  point  on  the  leading  edge  of  the 
profile.  The  stagnation  point  for  zero  incidence  is  located  a small 
distance  along  the  suction  side  from  the  tap  defined  as  the  leading  edge 
gauge(S./L  = 0.007).  This  is  why  the  Mach  number  at  S/L  = 0 does  not 
indicate  a stagnation  condition. 

On  the  pressure  side  the  velocity  distribution  for  all  exit  Mach 

numbers  is  similar.  There  is  a rapid  expansion  in  the  leading  edge  region. 

The  large  amount  of  turning  in  this  location  leads  to  an  over-expansion  of 

the  flow  around  S/L  «0.1.  After  a deceleration  in  this  region  there  is 
P P 

a favorable  pressure  gradient  to  the  trailing  edge  of  the  blade.  It  should 
be  noted  that  the  Mach  number  never  exceeds  unity  on  the  pressure  side  of 
the  blade.  This  is  expected  because  the  blade  design  is  such  that  the 
throat  occurs  at  the  trailing  edge  of  the  blade,  on  the  pressure  side  of 
the  profile. 

For  the  suction  side,  the  Mach  number  distributions  are  divided  into 
two  families  of  trends,  one  for  subsonic  exit  Mach  numbers  and  one  for  the 
low  supersonic  exit  Mach  numbers.  The  blade  was  designed  such  that  all  the 
flow  turning  would  be  accomplished  upstream  of  the  throat.  Demurentl] 
picked  this  scheme  in  designing  the  four  profiles  used.  The  designs  were 
such  that  the  effects  of  the  acceleration  due  to  flow  turning  would  be 
separated  from  the  acceleration  due  to  different  flow  areas.  In  the  region 
before  the  throat,  all  the  flow  turning  takes  place  and  after  the  throat 
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different  area  configurations  were  investigated.  It  is  the  flow  in  the 
region  after  the  throat  that  determined  the  two  families  of  trends  for  the 
reference  blade. 

For  the  subsonic  exit  Mach  numbers  the  flow  rapidly  accelerates 
from  the  stagnation  point  to  a peak  around  Ss/Lfi  “0.45.  This  over- 
expansion is  a direct  result  of  the  large  amount  of  turning  taking  place  in 
this  region.  Following  this  peak  the  flow  decelerates.  If  the  local 
velocity  is  supersonic,  this  velocity  reduction  occurs  through  a series  of 
Lambda  shocks  located  near  the  blade  surface.  For  exit  Mach  numbers  less 
than  0.7  the  velocities  are  subsonic  and  no  Lambda  shock  can  exist.  After 
this  region  the  flow  undergoes  a gradual  acceleration  to  the  trailing  edge 
for  the  case  of  subsonic  exit  Mach  numbers. 

In  the  transonic  regime,  the  flow  accelerates  from  the  stagnation 

point  to  a peak  at  S /L  *0.45.  Similar  to  the  high  subsonic  regime,  the 

s s 

isentropic  flow  velocity  is  reduced  by  a series  of  Lambda  shocks  located  on 

the  blade  surface  near  S /L_  s 0.5.  At  S /L_  »0.55  the  flow  decelerates 

s s s s 

to  unity,  to  meet  the  boundary  condition  of  choking  at  the  throat.  The 

flow  then  reaccelerates  rapidly  and  reaches  a second  peak.  At  this  point 

the  flow  passes  through  a shock  wave  and  the  isentropic  Mach  number 

decreases.  Inspection  of  the  Schlieren  photographs  shows  that  this  shock 

is  attached  to  the  trailing  edge  of  neighboring  blades.  The  intersection 

of  the  shock  and  the  suction  side  of  the  blade  occurs  further  downstream  as 

exit  Mach  number  increases.  This  is  the  reason  the  second  peak  in  Mach 

number  for  s 1.53  is  a little  further  along  the  suction  side  than  the 
ex 

peak  for  Mgx  s 1.09. 


-6- 


* 


W ' TT 

1.1.2  Positive  Incidence  a 

The  off-design  Mach  number  distribution  for  incidences  of  +1.5°  and 

♦15°  are  reported  in  Figures  3 and  respectively.  The  overall  trend  for 

the  positive  incidences  is  similar  to  that  for  zero  incidence,  but  some 

differences  do  exist.  On  the  pressure  side  the  flow  in  the  region  of  the 

leading  edge  still  has  a very  rapid  acceleration,  but  a smaller 

deceleration  occurs  in  the  region  of  S/L  -0.15  for  the  low  exit  Mach 

P P 

numbers.  As  the  exit  Mach  number  increases  an  overexpansion  does  take 
place  in  this  region.  This  overexpansion  gets  more  sevsre  as  exit  Mach 
number  increases,  but  disappears  at  the  highest  exit  Mach  number 
(M  *1.3).  For  both  positive  incidences,  at  M *1.3,  the  rapid 

6X  ©X 

acceleration  disappears  near  the  leading  edge  and  the  flow  undergoes  a 
gentle  acceleration  from  the  stagnation  point  to  the  trailing  edge.  One 
may  to  explain  this  is  that  due  to  the  high  curvature,  a separation  bubble 
is  attached  to  the  blade  surface  between  the  first  and  second  pressure  taps 
(Ss/Ls  a 0.07  and  Sa/La  * 0.17).  The  bubble  alters  the  streamlines  in  this 
region.  As  the  Reynolds  number  increases  (mass  flow  increases)  the 
boundary  layer  thickness  is  reduced  leading  to  a smaller  radius  of 
curvature  in  this  region.  This  increase  in  curvature  would  cause  the 
leading  edge  of  the  bubble  to  move  forward  until  at  the  highest  mass  flow 
it  is  over  the  first  pressure  tap.  The  forward  movement  of  the  bubble  is 
what  causes  the  sudden  change  in  the  shape  of  the  pressure  distribution 
near  the  leading  edge  as  mass  flow  increases. 

For  the  suction  side  no  noticeable  differences  exist  between  the  zero 
and  positive  incidence  data. 


1.1.3  Negative  Insldanaje 


For  the  negative  angle?  of  incidence,  significant  changes  occur  in  the 
Mach  number  distribution  (Figures  5 and  6).  The  stagnation  point  moves  up 
the  suction  surface  as  expected.  Due  tc  the  Large  amount  of  turning 
♦paired  to  accelerate  the  flow  from  the  stagnation  point,  around  the 
leading  edge  radius  to  the  pressure  side,  a rapid  increase  in  Mach  number 
i&  recorded  near  the  leading  edge  on  the  pressure  side  of  the  blade.  Again 
we  see  an  overexpansion  in  the  region  Sp/Lp  *0-1.  At  low  exit  Mach  numbers 
this  overexpansion  is  moderate,  followed  by  a smoothing  of  the  velocity 
distribifcion  in  the  mid-range  of  exit  Mach  numbers.  As  the  mass  flow  is 
increased  to  the  highest  mass  flow,  the  overexpansion  becomes  severe. 
Without  the  aid  of  flow  visualization  at  off-design  incidences  the 
following  explanation  cannot  be  verified,  but  heat  transfer  data  tends  to 
confirm  that  for  the  lower  mass  flows,  a long  separation  bubble  is  attached 
to  the  blade  surface  beginning  before  the  pressure  tap  at  Sp/Lp  s 0-17  and 
extending  over  the  next  three  pressure  taps.  This  is  indicated  by  these 
three  pressure  taps  reading  approximately  the  same.  At  the  higher  mass 
flows  the  separation  bubble  moves  forward  and  gets  smaller.  This  causes  a 
larger  overexpansion  followed  by  a more  rapid  reacceleration  after  the 
overexpansion.  The  bubble  collapses  due  to  the  acceleration  after  this 
point . 

On  the  suction  side,  there  is  a drastic  change  in  the  transonic  regime 
Mach  number  distribution  downstream  of  the  throat  for  the  negative 
incidences  when  compared  to  the  same  mass  flow  at  other  incidences.  For 
*11  mass  flows  the  region  from  th®  stagnation  point  up  to  the  throat  is 
similar  to  the  positive  and  zero  incidences,  and  for  the  subsonic  cases  the 

* 
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Mach  number  distribution  is  similar  from  the  throat  to  the  trailing  edge. 

For  supersonic  exit  conditions,  however,  the  flow  accelerates  further 

downstream  before  intersecting  with  the  trailing  edge  shock.  The  trailing 

edge  shock  intersection  with  the  suction  side  of  the  blade  can  be  seen 

progressing  downstream  as  the  exit  Mach  number  increases.  For  example,  the 

location  of  the  second  peak  in  isentropic  Mach  number  for  i * -15°  occurs 

at  S /L  « 0.70  for  M = 1.38,  but  for  M_„  t 1.70  this  peak  shifts  to 
s s ex  ex 

S /Ls  « 0.8.  The  same  is  true  for  the  -7.5°  incidence  except  that  at  the 
s 

highest  mass  flow,  the  limit  loading  condition  was  reached;  i.e.,  the 
trailing  edge  shock  did  not  intersect  the  suction  side  of  the  blade. 
Differences  were  also  caused  by  the  larger  inlet  area  of  the  negative 
incidence  cascade  geometry  which  changed  cascade  pressure  ratio. 

1.2  The  Effect  of  Incidence  on  Heat  Transfer 

The  calorimetric  heat  transfer  gauges  used  in  this  study  and  described 
in  previous  reports  were  designed  to  measure  heat  flux  as  the  hot  main  flow 
passes  over  the  blade  surface.  Heat  flux  cannot  be  measured  directly  but 
can  be  calculated  as  a function  of  the  change  in  temperature  of  the  gauge. 
The  short  run  time  of  the  experiment  allows  a guage  design  such  that  the 
heat  flux  is  proportional  to  the  temperature  rise  of  the  gauge  slug.  The 
change  in  temperature  of  the  slug  is  measured  with  a thermocouple.  The 
output  of  the  thermocouple  is  a constant  until  initiation  of  the  test,  then 
the  output  trace  starts  to  rise  montonically  until  the  test  is  terminated. 

The  gauge  is  modeled  as  a perfectly  insulated  cylindrical  disk  with 
one-dimensional  heat  transfer  into  the  slug  from  toe  hot  main  flow.  The 
heat  transfer  property  being  sought  is  the  convection  coefficient,  h, 
averaged  over  the  gauge  surface.  Since  the  heat  flux  into  the  gauge  can  be 
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Nisured,  the  oonveetion  coefficient  is  calculated  using  the  difference  in 
teaperature  of  the  hot  aain  flow  and  the  gauge  as  the  convective  driving 
force. 

1.2.1  Heat  Transfer  Distributions  Around  the  Blade 

The  heat  transfer  characteristics  of  the  reference  blade  were  plotted 
as  the  local  Nusselt  number  versus  the  blade  surface  location.  The  heat 
transfer  characteristics  were  measured  for  design  and  off-design  engine 
conditions.  The  results  are  generally  what  is  expected;  the  zero  and 
positive  Incidence  data  follow  the  same  trends  and  magnitudes,  while  the 
negative  incidence  data  has  slightly  different  trends  in  the  area  of  the 
loeading  edge  and  on  the  suction  surface  between  the  throat  and  the 
trailing  edge  (Figures  7 - 11).  These  two  regions  are  the  same  regions 
where  the  pressure  distribution  for  the  negative  incidence  varied  from  the 
zero  and  positive  incidences.  The  leading  edge  gauge  indicates  a lower 
value  of  heat  transfer  than  expected  in  the  region  of  the  stagnation  point. 
The  stagnation  point  is  not  centered  on  the  leading  edge  gauge  but  its 
location  is  within  the  gauge  area.  The  boundary  layer  grows  very  rapidly 
in  the  leading  edge  region  and  an  averaging  effect  is  expected  over  this 
region.  A dotted  line  is  shown  for  the  leading  edge  to  indicate  that  the 
actual  value  of  the  stagnation  point  is  expected  to  be  much  higher  than  the 
average  over  the  gauge  area.  The  trailing  edge  gauge  measures  heat 
transfer  from  all  sides  of  the  trailing  edge  radiius.  The  trailing  edge 
heat  transfer  rate  is  one  of  the  lowest  around  the  blade  surface.  The 
suction  side  laminar  boundary  layer  transition  zone  is  shown  at 
ii/L.  «0.2  by  a dotted  line.  The  nearest  gauges  were  centered  at 

9 3 


A 


-10- 


S_/L  • 0.13  and  0.24.  It  is  felt  the  transition  zone  must  occur  between 

these  two  beoause  other  likely  locations  had  a higher  density  of  gauges. 

For  zero  inoidence,  the  heat  transfer  to  the  pressure  side  is  exactly 

what  is  expected.  The  heat  transfer  decreases  rapidly  from  the  leading 

edge  as  the  boundary  layer  grows.  In  the  region  Sp/Lp  *0.2,  a slow 

transition  to  the  turbulent  boundary  layer  takes  place.  In  the  region 

after  the  transition,  the  heat  transfer  slowly  decreases  or  levels  off  to  a 

minimum  at  S /L  *0.5  and  increases  from  this  point  to  the  trailing  edge. 

P P 

Two  opposing  boundary  layer  properties  are  interacting  along  this  region. 

The  increased  length  would  indicate  that  the  boundary  layer  is  growing  in 
thickness,  while  acceleration  of  the  main  flow  would  cause  thinning  of  the 
boundary  layer.  An  inspection  of  the  velocity  distributions  indicate 
increased  acceleration  after  Sp/Lp  *0.5,  which  is  why  the  thinning  effect 
becomes  dominant  in  determining  the  heat  transfer  and  the  Nusselt  number 
rises  in  this  region. 

The  flow  on  the  suction  side  leads  to  a wide  variation  in  heat  j 

i 

transfer  as  a function  of  distance  along  the  suction  surface.  This  j 

variation  is  also  a strong  function  of  exit  Mach  number.  The  heat  transfer  j 

rate  decreases  from  the  stagnation  point  as  the  laminar  boundary  layer  j 

I 

grows.  The  laminar  transition  zone  occurs  in  the  vicinity  of  S /L  = 0.2.  j 

S 3 \ 

For  the  exit  Mach  numbers  of  0.62  and  0.79  there  is  an  increase  in  the  ' 

local  heat  transfer  rate  followed  by  a sharp  decrease  at  S_/L  = 0.38. 

S 3 

Inspection  of  the  pressure  distributions  indicate  that  this  is  where  the  1 

i 

peak  in  Mach  mmber  occurs.  In  the  following  region  of  rapid  flow  j 

deceleration,  a separation  bubble  forms  and  causes  the  reduction  in  heat 

transfer  to  the  blade.  After  this  the  heat  transfer  rate  tends  to  level 

off  for  the  remainder  of  the  blade.  For  the  exit  Mach  number  of  0.93,  this 

sharp  reduction  in  heat  transfer  occurs  at  S /L  *0.4  and  the  pressure  peak 

s s 
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occurs  at  S /L  -0.45.  There  appears  to  be  a "space  shift"  between  the 

heat  transfer  gauges  and  the  pressure  taps.  A possible  explanation  is  that 

the  heat  transfer  gauges  are  rough  enough  to  cause  the  transonic  flow  to 

undergo  deceleration  earlier  than  a smooth  surface.  This  is  suspected 

beoause  small  perturbations  of  surface  roughness  are  very  critical  in  the 

transonic  flow  regime.  At  S /L  -0.4  the  heat  transfer  rate  increasaes 

s s 

until  S./L.  *0.6.  This  corresponds  to  the  separation  bubble  reattaching 

9 8 

and  the  boundary  layer  thinning  after  reattachment.  After  this  point,  the 
heat  transfer  levels  off,  as  did  the  velocity  distribution. 

For  the  Hgh  mass  flows  (low  supersonic  exit  velocities)  the  suction 
side  heat  transfer  shows  a much  different  trend  downstream  of  the  throat. 
Again  minimum  heat  transfer  occurs  at  S /La  -0.40  and  the  decrease  in  the 
heat  transfer  rate  comes  before  the  velocity  peak.  This  systematic  "space 
shift"  between  the  velocity  maximum  and  the  heat  transfer  mininum  in  the 
transonic  regime  seems  to  indicate  that  gauge  roughness  is  a likely 
explanation.  After  this  minimum,  the  heat  transfer  increases  due  to 
reattachment  of  the  separation  bubble.  After  the  throat  (S  /L  ■ 0.55), 

9 9 

the  heat  transfer  continues  to  increase  with  acceleration  of  the  main  flow, 
which  thins  the  boundary  layer.  After  the  point  where  the  trailing  edge 
shock  system  intersects  the  blade  surface,  the  boundary  layer  separates. 

The  heat  transfer  decreases  and  levels  off  at  the  trailing  edge. 

For  the  positive  incidence  data  no  noticeable  difference  in  the  trend 
or  the  magnitude  of  the  heat  transfer  data  occurs.  The  effect  of  the 
instrumentation  "space  shift"  did  occur  and  indicates  the  heat  transfer 
gauges  do  alter  the  boundary  layer  in  the  transonic  flow  regime. 

For  the  negative  incidences  a few  differences  can  be  found  in  the  heat 
transfer  trends.  On  the  pressure  side,  the  slow  transition  of  the  laminar 
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boundary  layer  in  the  region  of  Sp/Lp  = 0.25  is  not  indicated  by  the  heat 
transfer  data  for  the  low  exit  Mach  numbers.  Instead  heat  transfer  is 
fairly  level  until  Sp/Lp  = 0.5,  where  rapid  acceleration  causes  thinning  of 
the  boundary  layer.  This  would  verify  the  long  separation  bubble  mentioned 
previously.  For  higher  mass  flows  the  heat  transfer  decreases  from  the 
leading  edge  gauge  and  then  increases.  At  the  highest  mass  flows  this 
decrease,  followed  by  an  increase  in  heat  transfer,  becomes  more  severe. 

As  discussed  previously,  as  the  mass  flow  increases,  the  separation  bubble 
becomes  smaller  in  length,  which  is  verified  by  the  decrease  and  downstream 
increase  in  the  heat  transfer  rate.  When  comparing  the  negative  incidence 
with  the  zero  and  positive  incidences  the  suction  side  heat  transfer  is 
very  similar  for  all  cases  except  the  low  supersonic  exit  velocities.  For 
these  high  mass  flow  runs  the  trailing  edge  shock  system  intersects  the 
suction  surface  further  down  the  blade  than  at  design  incidence.  After  the 
shock  hits  the  blade,  the  boundary  layer  separates  and  the  heat  transfer 
levels  off.  An  example  of  the  trailing  edge  shock  interaction  with  the 
boundary  layer  can  be  seen  by  comparing  the  exit  Mach  numbers  of  1.35  and 
1.71  for  the  negative  7.5°  incidence.  For  the  low  exit  Mach  number  the 
heat  transfer  levels  off  after  the  shock  system  intersects  the  suction 
surface,  but  for  the  high  exit  Mach  number,  where  limit  loading  was 
reached,  the  heat  transfer  continues  the  downward  trend  indicating 
thickening  of  the  boundary  layer. 


1.2.2 


The  Nusselt  nusber  averaged  over  the  blade  surface  was  calculated  for  1 

each- test  condition.  The  results  are  plotted  as  a function  of  exit  Mach 
number  for  different  incidences  (Figure  12).  The  results  shows  for  equal 

j ; 

aass  flow  rates  the  heat  transfer  for  off-design  conditions  is  within  10% 
of  that  for  design  incidence.  This  indicates  that  incidence  has  a very 
snail  effect  on  average  values  of  heat  transfer  for  the  conditions  tested. 

1.3  &— act  ftnd  Conclusions  on  the  Angle  of  Incidence  Studies 

The  angle  of  incidence  study  examined  the  effect  of  incidence  angle  on 
the  heat  transfer  performance  of  transonic  turbine  blades.  This  is 
iaportant  to  the  designer  who  must  be  assured  there  will  not  be  a 

. 

k catastrophic  loss  in  performance  at  off-design  angles  of  incidence.  It  is 

1 also  of  importance  to  the  turbine  designer  since  it  will  help  determine  the 

optimum  location  of  cooling  ports  in  order  to  avoid  overheating  of  the 

; \ 

blade  at  both  design  and  off-design  points. 

This  study  has  been  part  of  a larger  program  on  transonic  turbines 

which  has  been  going  on  at  the  M.I.T.  Gas  Turbine  Laboratory  for  the  past 

i 

j five  years.  The  aim  of  the  overall  program  is  to  gather  necessary  data  for 

\ cooled  transonic  turbine  blade  design  and  to  carry  out  the  testing  of  ; 

f i 

| film-cooled  transonic  blades. 
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1.3.1 


1.  The  off-design  incidence  angles  show  local  differences 
in  the  pressure  and  heat  transfer  distributions.  The 


I 

i 

i 


2. 


3. 


■ajor  areas  of  departure  from  the  design  case  are  on 
the  pressure  side  near  the  leading  edge  and  on  the 
suction  side  near  the  trailing  edge  - shock  interaction 

area. 

The  mean  Nusselt  number  was  found  to  be  independent  of 
the  incidence  angle  for  both  subsonic  and  low 
supersonic  exit  velocities.  For  equal  mass  flow  rates, 
mean  heat  transfer  rates  for  off-design  incidences  were 
within  10%  of  the  design  case. 

The  suction  side  laminar  transition  zone  was  found  to 
be  located  closer  to  the  leading  edge  than  reported  by 
Demuren[1]. 

The  pressure  distributions  reported  here  were  similar 
to  the  measurements  taken  by  Demuren  at  Von  Karman 
Institute  but  different  from  the  measurements  taken  at 
H.l.T.  The  heat  transfer  results  were  also  different, 
but  average  values  of  heat  transfer  were  in  good 
agreement.  Improvements  in  surface  instrunentation  and 
of  small  profile  changes  are  thought  to  be  the  cause  of 
the  difference  between  the  two  sets  of  data  taken  at 
the  M.l.T.  facility. 
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II.  Comparison  of  Test  Data  for  the  Reference  Convergent  Blade  with 
Straight  Bank  and  Thick  Trailing  Edge 
This  section  will  summarize  the  experimental  test  results  of  the  MIT 
wholly  convergent  turbine  blade  with  a straight  back  and  thick  trailing 
edge  conducted  in  transonic  cascade  tests.  Three  test  programs  of  interest 
will  be  summarized:  Harold  Demurer 's  work  at  VKI  and  MIT[1]  and  Fred 
Hajjar's  work  at  MIT[2],  In  the  MIT  studies,  both  blade  surface  static 
pressure  and  heat  transfer  distributions  were  measured,  whereas  the  VKI 
data  is  for  pressure  distribution  only. 

2.1  Siirnmarv  of  Demuren's  Work  at  VKI f 1 1 

2.1.1  The  VKI  Facility 

The  test  facility  used  for  the  turbine  cascade  investigation  was  the 
Von  Karman  Institute  high  speed  cascade  tunnel.  This  tunnel  is  of  the 
blowdown  type  and  is  supplied  with  high  pressure  air.  The  turbulence  level 
at  the  cascade  inlet  is  taken  to  be  0 . 7 Jfc . The  settling  chamber  pressure 
was  measured  with  a mercury  U-tube  manometer.  The  static  pressure  on  the 
cascade  blades  and  on  the  wind  tunnel  end  walls  were  measured  with  a 
mercury  multimanometer.  Mid-span  traverses  were  performed  at  the  inlet  and 
outlet  planes  of  the  cascade  utilizing  three  different  probes: 
at  the  inlet,  a pitot-directional  probe; 

at  the  outlet,  a pitot-directional  probe  combined  with  a single 
needle  static  probe  for  Mg  £.  1.3; 
at  the  outlet,  an  AVA-tube  probe  for  M^  > 1.2. 


2.1.2 


Schlieren  and  shadowgraph  photos  were  taken  at  various  exit  Mach 
ntsabers  to  verify  the  periodicity  of  the  outlet  flow  angle.  The  blade 
perfornanoe  was  determined  from  inlet  and  exit  traverses  and  blade  pressure 
distribution  measurements.  The  downstream  traverses  were  made  behind 
several  blades  including  the  instrumented  blades.  The  blade  pressure 
distributions  were  taken  without  probes  in  the  cascade  to  avoid 
probe-induced  disturbances . 

During  the  downstream  traverses  the  following  values  were  recorded 
continuously 

A PQ  * difference  between  the  settling  chamber  pressure 
and  the  total  pressure  at  the  probe. 

A P a difference  between  the  static  pressure  at  the  probe 
and  atmospheric  pressure. 

A PLR  a pressure  difference  measured  by  the  directional 
probe . 

The  inlet  Mach  number  was  based  on  the  total  pressure  and  area 
averaged  static  wall  pressure  measurements  0.1  chord  upstream  of  the 
cascade  inlet.  The  downstream  traverse  data,  taking  into  account  probe 
calibration,  was  used  to  calculate  an  efficiency  and  the  local  Mach  number 
Mg.  A calculated  Reynolds  number  was  based  on  the  inlet  flow  parameters 
and  referenced  to  1 cm  length. 


2.1.3  Blade  Geometry  and  lnatrumentation 


The  following  data  characterize  the  blade  and  cascade  geometry  for  the 
experiments  conducted  at  the  VK1  facility  (See  Figure  1): 

Blade  chord  c s 2.60  in. 

Blade  spacing  g = 1.95  in. 

Stagger  angle  = 51.0° 

Inlet  flow  angle  = 60.0° 

Trailing  edge  thickness  = .11  in. 

Number  of  blades  * 6 
g/c  * 0.75 

The  blade  pressure  distribution  was  measured  by  static  pressure  taps 
on  the  pressure  side  and  suction  side  of  two  neighboring  blades  (blades  3 
and  4)  such  that  the  instrumented  blade  surfaces  formed  a blade  flow 

passage. 

2.1.4  The  Measured. Elsa.  Field 

The  Mach  number  distribution  at  the  cascade  inlet  was  found  from  wall 
static  pressure  taps  in  a plane  .1  chord  ahead  of  the  leading  edge.  The 
inlet  Mach  number  variation  as  a function  of  exit  Mach  number  is  shown  in 
r'igure  13.  The  inlet  angle  variation  in  the  transverse  direction  was  also 
Tirasured  with  no  influence  of  the  exit  Mach  number  observed. 

The  blade  velocity  distribution  is  expressed  by  a local  Mach  number 
calculated  from  the  local  static  pressure  on  the  blade  surface  and  the 
total  pressure  upstream  of  the  cascade.  The  Mach  number  distribution  for  a 
geometry  g/c  = 0.75  is  shown  in  Figure  14. 
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2.1.5  The  Blade  Performance 

The  blade  performance  waa  evaluated  from  a downstream  wake  traverse  at 
an  axial  distance  of  .3  blade,  chord  behind  the  trailing  edge  of  the  blade. 
The  variation  in  blade  losses  as  a function  of  exit  Mach  number  is  plotted 


in  Figure  15.  For  exit  Mach  numbers  up  to  Mg  a 0.7  the  blade  losses  run  up 


to  8J  which  is  apparently  caused  by  the  trailing  edge  thickness.  A sudden 


rise  in  the  level  of  losses  from  M = 0.7  to  H : 0.85  can  be  attributed 

e e 


to  local  limbda  shock  losses  and  a shock-boundary  layer  separation  along 
the  upper  suction  side  of  the  blade.  The  decrease  in  losses  from  11%  at 


M a 0.85  to  around  5)1  at  M.  s 1.3  (the  design  point)  can  be  attributed  to 


the  shocks  becoming  oblique  and  flow  reattachment  occurring.  It  is 


suggested  that  the  increase  in  losses  for  Mg  > Mdesign  is  due  to  the 


increasing  strength  of  the  adjacent  blade  left  running  trailing  edge  shock 
which  intersects  the  suction  side  and  causes  the  boundary  layer  to 
separate.  The  losses  recorded  for  Mfi  > load  depend  on  the  total 

pressure  losses  due  to  the  boundary  layer,  the  trailing  edge  shocks,  and 
the  mixing  processes  between  the  exit  and  the  traversing  planes. 

The  outlet  flow  angle  8g  as  a function  of  exit  Mach  number  is  also 


shown  in  Figure  15.  the  flow  angle  decreases  slightly  from  25.5°  to  24° 


betweem  Mft  = 0.6  and  1.3.  For  Mg  > Mdesign  design  = 1>3)’  the  flow  exit 


angle  increases  slowly  up  to  Mg  = 1.4,  but  then  increases  sharply  for 


Me  > 1.4.  This  sharp  increase  can  be  associated  with  an  over-expansion  at 


the  blade  trailing  edge.  At  Mg  = load  (M  limit  load  = 1.59)  the 


outlet  flow  angle  deviation  was  about  10 


1 
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2.1.6 


Tests  were  also  conducted  to  determine  the  effects  of  stage  solidity 
(chord  length/blade  spacing)  on  the  blade  performance.  The  blade  was 
tested  at  g/c  = 0.81  and  0.695.  Examination  of  the  resultant  blade 
velocity  distributions  and  the  SchHeren  photos  showed  a change  in  the 
locations  and  inclinations  of  the  left-running  shocks  for  comparable  exit 
Mach  n embers.  Also,  as  the  blade  spacing  increased,  the  shock  strength 
became  weaker  and  so  flow  separations  resulting  in  the  shock-boundary  layer 
interactions  were  avoided.  No  other  significant  effect  on  the  blade 
velocity  distribution  was  seen. 


i 

j 


2.1.7  Effect  of  Reynold's  Number  Variation 

The  results  of  a series  of  tests  to  examine  the  effects  of  Reynold's 

5 6 

number  variation  between  10^  and  10  showed  no  noticeable  change  in  either 

the  blade  pressure  distribution  nor  in  the  shock  structure.  Hence  it  was 

5 6 

concluded  that  Reynold's  number  effects  in  the  range  of  10^  and  10  showed 

no  noticeable  change  in  either  the  blade  pressure  distribution  nor  in  the 

shock  structure.  Hence  it  was  concluded  that  Reynold's  number  effects  in 
5 6 

the  range  of  10  to  10  were  negligible. 


r 
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• i<8  Earn  .Stream.  Mato 

Analysis  of  the  downstream  wake  using  the  downstream  traverse  .3  chord 
behind  the  blades  shows  an  interesting  effect  of  Mach  number  on  the. 
structure  of  the  wake.  The  downstream  wake  profile  is  shown  in  Figure  16. 
As  the  exit  Mach  number  increases,  the  flow  asymmetry  increases.  A 
characteristic  effect  is  the  increase  in  the  wake  depth  during  transition 
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from  subsonic  Mach  number  to  supersonic  Mach  number.  The  width  of  the  wake 
also  changes.  Careful  examination  of  the  Schlieren  photographs  of  the  wake 
flow  shows  that  it  consists,  under  certain  conditions,  of  von  Karman  vortex 
sheets . 

2.2  Siaamrv  of -K.I.T.  Work 

Deauren's  and  Hajjar's  work  was  conducted  at  the  M.I.T.  Gas  Turbine 
Laboratory  utilizing  the  hot  blowdown  turbine  facility.  This  facility 
delivers  air  at  selected  pressures  and  temperatures  for  a maximum  duration 
of  ten  seconds  at  maximum  flow  rate.  Tests  have  shown  that  steady-state 
flow  conditions  are  reached  0.2  seconds  into  the  test. 

2.2.1  Deauren's  M.I.T.  Workfll 

Demuren's  experiment  at  M.I.T.  tried  to  duplicate  most  of  the  test 
parameters  utilized  in  the  VKI  tests.  The  primary  differences  between  the 
two  were: 

1.  g/c  a 0.695  whereas  VKI  g/c  = 0.75 

2.  9 blade  cascade  whereas  VKI  utilized  6 blades 

3.  short  duration  whereas  VKI  ran  continuously  steady  state 

4.  10H  turbulence  level  vs.  0.7%  at  VKI 

5.  flow  temperature  of  400°  F vs.  room  temperature  at  VKI. 

The  same  blade  profile  was  utilized  in  both  tests.  Also,  the  blade 
was  instrumented  in  the  same  manner. 
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2.2. 1.1 


The  plot  of  the  measured  blade  velocity  distribution  is  shown  in 
Figure  17.  A very  fast  expansion  is  seen  to  take  place  on  the  suction  side 
with  the  Mach  number  reaching  a peak  at  a location  where  some  surface  bumps 
are  located  on  the  blade,  the  flow  decelerates  a bit  and  then  quickly 
starts  to  accelerate  again  until  the  surface  is  intercepted  by  the  left- 
running trailing-edge  shock  from  the  adjacent  blade,  downstream  of  which 
the  flow  decelerates  and  later  accelerates  towards  the  trailing-edge.  As 
the  pressure  ratio  increases  across  the  cascade,  the  intercepting  shock 
moves  down  towards  the  blade  trailing  edge.  The  limit  loading  condition  is 
such  that  this  shock  passes  downstream  of  the  trailing  edge  and  no  longer 
affects  the  channel  flow.  The  flow  on  the  pressure  side  rapidly 
accelerates  immediately  downstream  of  the  stagnation  point,  quickly 
decelerates,  and  then  gradually  accelerates  to  its  maximum  value  at  the 
trailing  edge.  Once  the  flow  is  choked,  there  is  very  little  change  in  the 
pressure  side  velocity  distribution. 

2. 2. 1.2  Comparison  Between  Mil  and  VKI  Pata 

A detailed  comparison  was  madetl]  cf  the  Mach  nunber  distribution 
around  the  blade  obtained  at  the  MIT  hct  blowdown  facility  with  that 
obtained  at  the  VKI  facility.  Overall,  the  pattern  of  the  Mach  number 
distribution  was  similar.  Fast  acceleration  from  the  stagnation  point  on 
| the  suction  side  led  to  the  first  Mach  number  peak.  But  the  location  and 

s 

? absolute  magnitude  of  the  peak  were  different.  The  peak  was  shifted 

forward  in  the  MIT  data.  Several  possible  reasons  were  suggested  for  this 
shift.  The  VKI  tests  were  conducted  with  cold  flow  and  a low  turbulence 
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level  (.7$)  whereas  the  MIT  work  was  with  hot  flow  air  (400°  F)  with  a 
turbulenoe  level  of  101.  Schlieren  photos  taken  at  VKI  show  a distinct 
separation  region  on  the  suction  side  of  the  blade.  The  combined  effects 
of  large  blade  eurvature  and  surface  roughness  could  have  triggered  the 
separation.  The  high  free-stream  turbulence  in  the  MIT  facility  could 
allow  a more  rapid  flow  reattachment  and  thus  change  the  velocity  profile 
in  this  area.  Also  minor  changes  in  static  pressure  tap  locations  on  the 
blade  surfaces  could  account  for  shifts  in  the  peak  locations  and 
magnitudes.  The  overall  pattern  of  the  velocity  distribution  is  the  same 
in  both  sets  of  data.  Good  agreement  is  seen  on  the  pressure  side 
distribution  with  a gradual  Mach  number  increase  up  to  the  trailing-edge. 


I 

| 
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2.2. 1.3.  Dawuren's  Blade  Heat  Transfer  Distributions 

The  measured  blade  surface  heat  transfer  distributions  are  presented 
in  Figure  18.  Overall,  the  blade  surface  velocity  and  cascade  turbulence 
level  greatly  Influence  the  boundary  layer  properties  which,  in  turn, 
influence  the  heat  transfer  rate  to  the  cascade  blade  surface.  A similar 
effect  on  both  the  velocity  and  heat  transfer  distributions  can  be  seen  in 
the  area  of  the  shock-boundary  layer  interaction.  The  rapid  acceleration 
from  the  leading  edge  on  the  suction  side  creates  a condition  very 
favorable  for  the  formation  of  a laminar  boundary  layer  and  a resultant 
decrease  in  Nusselt  number.  The  transition  to  turbulent  boundary  layer 
flow  causes  a sudden  increase  in  the  Nusselt  number.  Thereafter  the 
Nusselt  number  fluctuates  in  a manner  similar  to  the  velocity  distribution 
until  it  intercepts  the  left-running  shock  from  the  adjacent  blade. 
Downstream  of  the  shock,  the  Nusselt  number  decreases  significantly  and 
then  increases  as  the  flow  accelerates  toward  the  trailing  edge. 
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On  the  pressure  side  of  the  blade,  the  rapid  acceleration  from  the 
stagnation  point  is  again  conducive  to  the  formation  of  a laminar  boundary 
layer  with  the  resultant  decrease  in  Nusselt  number.  Transition  appears  to 
occur  at  x/c  > 0.3  with  the  Nusselt  number  increasing  to  the  trailing  edge 
due  to  the  continuous  acceleration  of  the  flow.  This  increase  in  the 
Nusselt  number  is  attributed  to  the  thinning  of  the  boundary  layer  as  the 
flow  accelerates.  A high  level  of  heat  transfer  was  recorded  close  to  the 
trailing  edge  (about  75S  of  the  mean  heat  transfer  measured  at  the  leading- 
edge  area).  Because  of  the  finite  size  of  the  heat  transfer  gauge 
(3/32  in.  diameter),  the  heat  transfer  rate  is  averaged  over  this  area, 
whereas  the  Nusselt  number  could  have  locally  varied  widely  in  this  area. 
Thus  the  local  heat  transfer  rates  are,  in  fact,  average  heat  transfer 
rates  over  a given  small  area. 

2.2.2  Hfl.Uftr'fl  MIT  Hart[2] 

Hajjar' a work  was  intended  to  first  correlate  the  work  of  Demuren  and 
then  investigate  the  angle  of  incidence  effects  on  the  cascade  blade 
performance.  The  same  blade  profile  was  utilized  except  Hajjar  removed  the 
surface  bumps  referred  to  by  Demuren  in  his  thesis.  The  same  principles  of 
instrumentation  were  used  by  both  Hajjar  and  Demuren.  Hajjar' s results  are 
presented  in  Section  1 of  this  report. 

; 

2.3  Comparison  of  Data 

The  Mach  number  distribution  measured  by  Hajjar  around  the  reference 
blade  at  zero  angle  of  incidence  (Figure  19)  was  compared  to  the  data 
published  by  Demuren[1],  Demuren  presents  two  different  sets  of  Mach 
number  distributions,  one  for  data  recorded  at  M.I.T.  utilizing  a hot  flow 
| (Figure  17)  and  one  for  data  recorded  at  V.K.I.  utilizing  a room 
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temperature  flow  (Figure  14),  for  basically  the  same  reference  blade. 

Hajjar’s  and  Demuren's  N.l.T.  data  shows  good  agreement  on  the 
pressure  aide  of  the  profile  and  poor  agreement  on  the  suction  side.  On 
the  suction  aide,  Deauren'a  first  Mach  number  peak  occurs  for  a smaller  x/c 
and  ia  lower  in  magnitude.  Demuren's  second  peak  for  the  low  supersonic 
exit  Mach  number  ia  also  shifted.  Comparing  Demuren's  V.K.l.  data  with 
HaJJar's  M.I.T.  testa  shows  much  better  agreement.  Demuren  discussed 
reasons  such  as  the  suotion  side  surface  profile  "bump",  turbulence  level 
(M.I.T.  10K  vs.  V.K.l.  0.7%)  and  the  difference  in  the  mainstream  flow 
temperature  to  account  for  the  difference  in  his  two  sets  of  data.  With 
instrumentation  improvements  discussed  by  Ha j jar [2],  the  V.K.l.  data  was 
reproduced  at  the  M.I.T.  facility.  It  is  possible  that  the  roughness 
oaused  by  instrumentation  on  Demuren's  M.I.T.  blade  led  to  deceleration  of 
the  transonic  flow  prematurely.  This  would  account  for  the  velocity  peaks 
occurlng  sooner  in  Demuren's  data. 

A sample  comparison  of  heat  transfer  was  made  with  Demuren's  work 
(Figure  20).  The  teat  conditions  compared  were  an  exit  Mach  number  of  1.09 
for  Ha J Jar' a tests  to  an  exit  Mach  number  of  1.08  for  Demuren's  data. 

There  are  differences  in  the  local  Nusselt  number  in  all  sections  of  the 
blade  except  the  fully  turbulent  region  after  the  throat  on  the  suction 
side.  Since  the  pressure  distributions  were  much  different  in  this  region, 
one  should  not  expect  a good  correlation  of  heat  transfer  results.  For  the 
pressure  side  much  better  agreement  was  expected,  but  the  results  were 
somewhat  different.  An  explanation  is  that  Demuren's  gauges  were  in  one 
axial  line,  which  may  cause  boundary  layer  perturbations  as  discussed  by 
Ha J Jar. 
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Plows  in  gas  turbines  are,  in  general,  unsteady.  The  unsteady  effects 
arise  from  factors  such  as  blade  passage,  propagation  of  azimuthal 
nonuniforaities , and  the  propagation  of  periodic  disturbances  which 
originate  upstream  of  the  turbine.  These  unsteady  disturbances  are  likely 
to  have  a narked  effect  on  the  blade  passage  flow  field,  the  base  flow,  the 
trailing  edge  shock  system,  and  the  blade  heat  transfer.  The  aerodynamic 
flow  around  the  trailing  edge  is  of  great  interest  because  it  strongly 
influences  blade  losses. 

When  a transonic  cascade  is  operating  at  conditions  other  than  limit 
loading,  the  left  running  trailing  edge  shock  of  one  blade  usually 
interacts  with  the  boundary  layer  of  the  suction  side  of  the  adjacent 
blade,  and  this  interaction  can  result  in  local  or  full  flow  separation, 
which  will  have  an  influence  on  the  turbine  flow  losses.  Since  the  rotor 

typically  operates  at  a high  angular  velocity  with  many  individual  blades, 

5 

a stator  can  see  in  the  neighborhood  of  10  fluctuations  per  minute.  The 
purpose  of  this  investigation  is  to  examine  the  effects  of  high  frequency 
periodic  unsteadiness  on  the  performance  of  a given  turbine  blade  by 
simulating  the  blade  passage  disturbance  in  a 2-D  transonic  turbine 
cascade. 

The  investigation  will  examine  the  following  phenomena: 

1 . The  aerodynamic  and  heat  transfer  properties  in  the 
neighborhood  of  the  oscillating  shock-boundary  layer 
interaction  on  the  suction  side  of  an  adjacent  blade; 

2.  The  effect  of  periodic  unsteady  disturbances  on  the 
potential  flow  in  the  blade  passage; 

3.  The  effect  of  oscillating  disturbances  on  the  base 
pressure  behind  the  stationary  blade; 

4.  The  accuaulated  losses  in  the  wake 
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The  investigation  will  help  provide  a basic  understanding  of  the  loss 
aechaniaas  which  are  inherent  in  all  gas  turbine  generators  which  utilize 


transonic  rotating  machinery. 

The  experiment  will  be  performed  in  the  M.I.T.  linear  cascade  tunnel 
attached  to  the  hot  blowdown  facility.  A high  frequency  disturbance 
generated  by  rotating  a small  elliptical  body  downstream  of  the  turbine 
cascade  row  will  be  used  to  simulate  the  periodic  unsteady  effects.  Blade 
surface  and  cascade  wall  instrumentation  will  be  utilized  in  measuring  the 
aerodynamic  and  heat  transfer  properties  for  different  operating 
conditions.  A traversing  probe  will  be  used  downstream  to  measure  the 
cascade  wake  flow  properties.  It  is  anticipated  that  a model  will  be 
developed  as  a result  of  the  experimental  investigation  which  will  include 
the  effects  of  the  shock-boundary  layer  interaction,  the  base  pressure 
fluctuations,  and  the  resultant  losses  in  the  wake. 

To  date,  the  air  turbine  used  to  rotate  the  elliptical  rotor  has  been 
extensively  tested,  the  cascade  has  been  modified  to  accept  the  unsteady 
experiment  and  the  pressure  field  around  the  rotor  has  been  mapped. 

3.1  The  -Elowfleld  Generated  bv  a Rotating  Ellipse 

3.1.1  Assumptions 

For  this  study  a two-dimensional,  inviscid,  incompressible, 
lrrotatlonal  fluid  was  assumed.  The  rotating  elliptical  body  is  an  ellipse 
described  by  the  equation: 

2 2 

-V-c.i 

a2  b2 

where  a and  b are  the  semi-major  and  semi-minor  axes  respectively. 
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Ih«  complex  potential  which  describes  the  fluid  motion  generated  by  a 
rotating  elliptical  cylinder  about  its  axis  in  an  infinite  fluid  is  given 
in  terms  of  elliptical  coordinates  as 

W s i C e"2J  (2) 

with  the  velocity  potential  given  as 

* s Ce_2e  sin  2r> 

and  the  stream  function  given  as 

t s Ce“2t  cos  2n 

with  C given  as 

C = 1/M  u(a  + b)2  (5) 


where  w is  the  angular  velocity  of  the  rotating  ellipse  and  e and  n are 
elliptical  coordinates. 

3.1.3  Elliptical  Coordinates  

The  elliptical  coordinates  (e,  n.  z)  form  an  orthogonal  curvilinear 
coordinate  system  which  is  related  to  a cartesian  coordinate  system  (x,y,z) 


x s c cosh  t cos  n 
y = c sinh  sin  n 


where  0 i £ <00  and  0 £ n <2*.  curves  of  constant  e give 

. f . , 

c*  cosh2e  c^sinh^e 


which  are  con focal  ellipses  with  foci  at  (±  c,0).  Curves  of  constant  n 


\ 

f. 

I 

j 

i 

\ 

1 


give 


2 2 2 2 
c cos  n c sin  n 


= 1 


(10) 


which  are  confocal  hyperbolas  with  foci  at  (±  e, 0). 


3.1.11  Flow  Field  Velocity  Components 

The  velocity  vector  J L can  be  resolved  into  elliptical  coordinate 
components 

1 = V u + V_  u_  (ID 

e e m n 

where  uc  and  un  are  unit  vectors.  The  velocity  components  are  given  by 


Ve 


. _1 a_ 


he  3 e 


and 


_ .Ml 
n hn 


v-  " K 3 n 


with  the  scale  factors 


12  2 

h a h = c Vcosh  e - cos  n 
e n 

substituting  equation  (3)  and  (14)  into  (12)  gives 


-2  <t> 

V ■ /'  1 " 2 2 

e -yc  cosh  e - cos  n 


and  equations  (3),  (4),  and  (14)  into  (13)  gives 


2i_ 


V - / “2 T 

n yc  cosh  e - cos  r 


(12) 


(13) 


(14) 


(15) 


(16) 
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Following  the  analysis  by  Milne-Thomson,  an  expression  for  the 
pressure  and  velocity  in  the  slipstream  of  an  advancing  propeller  can  be 


written  as: 


£ ♦ 1/2  V2  - r2  u fi 


where  P is  the  local  pressure,  p the  density,  V the  absolute  air  speed, 
in  the  angular  velocity  of  the  rotating  body,  n the  angular  speed  of  a plane 
containing  the  point  of  interest,  r the  distance  containing  the  point  of 
interest,  and  « the  pressure  at  infinity.  Using  equation  17,  the  pressure 
field  can  be  determined  using  the  results  of  the  potential  velocity  field. 


In  order  to  better  understand  the  aerodynamics  of  the  elliptical 
rotor,  a test  has  been  conducted  to  measure  the  flowfield  generated  by  the 
rotating  body.  Three  small  (.24  in.  diameter  diaphragm)  miniature  quartz 
Kistler  Series  600  pressure  transducers  were  placed  with  the  pressure 
sensitive  diaphragms  parallel  and  flush  to  a coverplate  to  record  the  short 
term  static  pressure  responses  generated.  Figure  21  shows  the  location  of 
the  transducers  in  relation  to  the  rotating  body.  Transducer  location  1 
corresponds  to  the  mean  radius  of  the  elliptical  body.  Location  2 
corresponds  to  one  mean  radius  past  the  maximum  radius  of  the  rotor, 
whereas  location  3 is  approximately  2 1/4  mean  radii  past  the  maximum 
radius. 

Figure  21  details  the  specifications  of  the  elliptical  body  where  the 
semi  major  axis  a is  0.632  inches,  the  semi  minor  axis  b is  0.211  inches, 


r 


and  the  focii  are  at  c * ± 0.596  inches.  In  the  analysis  it  is  assumed 
that  this  body  is  described  by  equation  (1),  an  ideal  ellipse.  Utilizing 
the  model  just  described,  flow  field  calculations  were  made  using  the  three 
transducer  locations  as  points  of  interest.  The  first  location,  however, 
at  r « 0.42  inches  was  moved  to  r s 0.632  inches  in  the  computations  to 
avoid  singularities  in  the  calculation. 

The  position  of  the  elliptical  body  relative  to  the  three  points  of 

interest  was  described  by  a 24  "quasi-steady"  point  coordinate  system  at 

rest  with  respect  to  the  rotating  body.  All  is  symetric  about  180°.  See 

Figure  22  and  Table  2.  In  order  to  determine  the  velocity  potential 

(equation  3)  at  each  point  of  interest  the  elliptical  coordinates  for  each 

point  were  determined  using  equations  6,  7,  and  9 in  an  iterative 

procedure.  The  computed  elliptical  coordinates  are  tabulated  in  Table  2. 

Equations  3i  4,  and  5 are  used  to  determine  the  value  of  potential  and 

stream  function  for  each  point  of  interest.  The  constant  C (equation  5) 

2 

was  determined  to  be  6.5  ft  /sec  for  angular  rotation  of  50K  RPM  and 

2 

10.3  ft  /sec  for  80K  RPM.  Representative  values  of  $ and  Y are  presented 
in  Table  3 and  plotted  in  Figures  23  and  24.  the  resultant  velocities  were 
calculated  using  equations  15  and  16  and  are  presented  in  Table  4 and 
Figures  25,  26  and  27. 

Equation  17  was  used  to  calculate  the  flowfield  pressure  at  the 
various  points  of  interest.  Since  in  this  study  we  are  interested  only  in 
changes  of  pressure  for  a given  rotation  of  the  elliptical  body,  the 
equation  can  be  simplified  to 

P1  - Pj  = 1/2  p (Vj2  - VA2) 
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Table  5 presents  the  results  of  the  pressure  calculations  in  terms  of 

dynamic  pressure.  As  can  be  seen,  the  pressure  change  at  location  1 is  'i . 

i 

i . 

much  greater  than  at  locations  2 and  3 for  both  cases  under  consideration. 

In  terms  of  pressure  changes  per  revolution,  at  location  1 the  total  change 

is  1.10  psi  at  50K  RPM  whereas  at  80K  RPM  the  change  is  2.81  psi.  The  j 


following  summarizes 

the  calculation: 

Location 

50K  RPM 

80K  RPM 

1 

1.10 

psi 

2.81 

psi 

2 

1.55 

x 10~3 

psi 

3.98 

x 10"3 

psi 

3 

.88 

X 

o 

1 

VJ1 

psi 

2.24 

x IQ"5 

psi 

The  results  of  the  air  turbine  tests  for  transducer  location  1 (under 
the  elliptical  body)  are  shown  in  Figure  28  and  in  Table  6.  In  the  graph 
1 psi  is  about  27  mv.  Hence,  for  50K  RPM,  the  average  output  was  .47  psi 
whereas  for  80K  RPM,  1.1  psi.  The  analogous  procedures  for  calculating 
pressures  for  20K,  30K,  40K,  60K,  70K,  and  90K  RPM  have  been  performed  and 
the  results  are  presented  in  Figure  29.  Here  the  results  of  calculations 
and  measurements  are  presented.  The  maximum  aP  values  for  each  RPM  are 
recorded,  and  show  a closer  correlation  to  the  theoretical  results. 
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DESIGNATION 

X-COORD 

Y-COORD 

S/L  ! 

• ; 

Pressure  Side  and  Leading  Edge 

! 

PI 

2.415 

0.055 

0.92 

P2 

2.300 

0.110 

0.87 

P3 

2.180 

0.155 

0.82 

P4 

2.015 

0.210 

0.76 

P5 

1.815 

0.265 

0.68 

P6 

1.630 

0.305 

0.59 

P7 

1.315 

0.345 

0.49 

P8 

1.030 

0.340 

0.38 

P9 

0.770 

0.300 

0.28 

P10 

0.500 

0.195 

0.17 

Pll 

0.285 

0.040 

0.07 

LE 

0.115 

0.010 

0 f 

Suction  Side  and  Trailing  Edge 

1 

i 

i 

S17 

0.005 

0.165 

0.06 

S16 

0.005 

0.410 

0.13 

S15 

0.115 

0.705 

0.22 

S14 

0.230 

0.830 

0.26 

S13 

0.330 

0.900 

0.30  1 

S12 

0.440 

0.945 

0.33 

Sll 

0.550 

0.965 

0.36 

S10 

0.740 

0.945 

0.41 

S9 

0.830 

0.920 

0.44 

S8 

1.070 

0.820 

0.51 

S7  .. 

1.290 

0.720 

0.58  j 

S6 

1.535 

0.610 

0.65 

S5 

1.745 

0.510 

0.72  : 1 

S4 

1.970 

0.405 

0.79  ! i 

S3 

2.185 

0.300 

0.85 

S2 

2.400 

0.195 

0.92 

SI 

2.535 

0.120 

0.96 

TE 

2.605 

0.030 

1 ■ ! 

Dimensions  given  in  inches,  T ■ *0.005 


TABLE  1.  LOCATION  OF  PRESSURE  TAPS 
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FIGURE  so.  LOCAL  NUSSELT  NUMBER  VS  BLADE  SURFACE  LOCATION 


FIGURE  fttf  LOCAL  NUSSELT  NUMBER  VS  BLADE  SURFACE  LOCATION 
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FIGURE  10R  LOCAL  NUSSELT  NUMBER  VS  BLADE  SURFACE  LOCATION 
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FIGURE  LOCAL  HUSSELT  NUMBER  VS  BLADE  SURFACE  LOCATION 


of  Incidence  • -7.5 


FIGURE  10D  LOCAL  NUSSELT  NUMBER  VS  BLADE  SURFACE  LOCATION 
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FIGURE  11A  LOCAL  NUSSELT  NUMBER  VS  BLADE  SURFACE  LOCATION 


FIGURE  JIJC  LOCAL  NUSSELT  NUMBER  VS  BLADE  SURFACE  LOCATION 
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FIGURE  20  COMPARISON  OF  HEAT  TRANSFER  DATA 
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Figure  24  Stream  Function  for  Transducer  Location 


Figure  25  Velocity  Calculation  at  Transducer  Location 
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Figure  27  Velocity  Calculation  at  Transducer  Location  3 
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